Engineering the thermopower in semiconductor-molecule junctions: towards high 

thermoelectric efficiency at the nanoscale 



D. Nozaki, H. Sevincli, W. Li, R. Gutierrez, and G. Cuniberti 

Institute for Materials Science and Max Bergmann Center of Biomaterials, 
Dresden University of Technology, 01062 Dresden, Germany 
(Dated: June, 05 2009) 

We propose a possible route to achieve high thermoelectric efficiency in molecular junctions by 
combining a local chemical tuning of the molecular electronic states with the use of semiconducting 
electrodes. The former allows to control the position of the HOMO transmission resonance with 
respect to the Fermi energy while the latter fulfills a twofold purpose: the suppression of electron- 
like contributions to the thermopower and the cut-off of the HOMO transmission tails into the 
semiconductor band gap. As a result a large thermopower can be obtained. Our results strongly 
suggest that large figures of merit in such molecular junctions can be achieved. 
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Introduction— Thermoelectric materials convert ther- 
mal gradients and electric fields for power conversion and 
for refrigeration, respectively. With increasing energy de- 
mand, thermoelectric applications are attracting a con- 
siderable interest. Unfortunately, thermoelectrics find 
currently only special applications due to their limited 
efficiency, which is measured by a dimensionless param- 
eter, the thermoelectric figure of merit: ZT = S'^TG/tv, 
which includes the Seebeck coefficient (thermopower) 
an average temperature T, the electrical conductance G, 
and the thermal conductance = A^ei + /^ph, the latter 
containing both electronic A^ei and vibrational A^ph con- 
tributions. Maximizing ZT is challenging, because op- 
timizing one physical parameter often adversely affects 
another. 

However, as suggested in the early 90s [^], dimension- 
ality reduction towards the nanoscale could provide an 
additional parameter to tune the electrical and ther- 
mal response of thermoelectric materials. This has 
triggered an active research on new nanoscaled ther- 
moelectric materials [2, 3, 4, 5, 6, 7, 8, 9, 10, 11]. 
An alternative to inorganic-based materials could be 
to exploit molecules, which have been already exten- 
sively investigated in the context of charge transport 
and molecular electronic applications in the past two 
decades [ ]. Indeed, the Seebeck coefficient of differ- 
ent single molecules [ ] as well as the thermal conduc- 
tance of self-assembled monolayers [ ] have been mean- 
while experimentally investigated. Theoretically, only 
few studies on the thermoelectric properties in few level 
systems (molecules, small quantum dots) have been pre- 
sented to date [15, 16, 17, 18, 19, 20, 21, 22, 23, 24]. Con- 
cerning this "organic route" , two aspects should be taken 
into account. First, as proposed by Mahan and Sofo [ ], 
the presence of a sharp resonance near the Fermi level 
Ep can considerably increase the thermopower, since the 
latter depends on the derivative of the conductance near 
Ep. One key advantage of using molecules as potential 
thermoelectrics is the capability to tune their chemical. 



and hence also their electrical and thermal properties in 
a very controlled way. The second aspect is more generic: 
by using metallic electrodes to contact a molecule, there 
are hole-like and electron-like contributions to S aris- 
ing from tunneling through the highest-occupied molec- 
ular orbital (HOMO) and lowest-unoccupied molecular 
orbital (LUMO) states, respectively. Both contributions 
appear with different signs in S and thus partially cancel 
each other. Could one overcome this problem? 
In this Letter, we propose a combination of chemical tun- 
ing of the molecular electronic structure with the use of 
semiconducting electrodes to improve the thermoelectric 
efficiency of molecular junctions. Such electrodes pro- 
vide a way to lift the near cancelation of hole-like and 
electron-like contributions to the thermopower. Specifi- 
cally, if the LUMO resonance lies within the semiconduc- 
tor band gap, no states are available for the "electronic" 
channel, which is then blocked and the thermopower 
will be mainly determined by the HOMO (hole) chan- 
nel. This together with the inclusion of donor groups in 
a TT-conjugated molecule (chemical tuning) can lead to 
a dramatic increase of the thermopower of the junction. 
We illustrate these ideas by using a first-principle based 




FIG. 1: A typical molecular junction investigated in this 
study. Two silicon electrodes terminated at the (111) facet are 
bridged by a molecule. In the calculations, periodic bound- 
ary conditions in the lateral directions have been used. The 
electrode surface is passivated with hydrogen atoms to avoid 
strong structural distortions related to surface reconstruction 
effects. 
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approach to the electronic structure and the transport 
properties of the molecular junction as well as by a min- 
imal model Hamiltonian. A typical electrode-molecule 
set-up used in our simulations is displayed in Fig. 1 [26]. 
Effective model Hamiltonian— We will first illustrate 
our approach within a minimal model Hamiltonian in 
order to highlight different factors influencing the ther- 
mopower of the junction. We consider two electronic 
levels, which mimic the frontier orbit als of a molecule 
and which are coupled to semiconducting electrodes as 
shown in Fig. 2(a). The thermopower S and the fig- 
ure of merit ZT [2/] are calculated as a function of the 
relative position of the HOMO to the Fermi level Ef^ 
A = £^HOMO — Ef^ and for different coupling strengths 
to the electrodes. In Fig. 2(b) and (c), room temperature 
values for S and ZT are shown. We see that it is pos- 
sible to obtain S values as large as ^2 mV/K when the 
molecule-electrode coupling is ~10 meV, and the HOMO 
level is placed 0.5 eV below Ep- Even much higher val- 
ues of S are possible for weaker coupling strengths. The 
optimum value of A is related to the operating tem- 
perature and it is smaller for lower temperatures [ ]. 
Very high thermopower can be achieved with a delta- 
function shaped conductance peak near to Ef for bulk 
systems [ ] ; for a model molecular system this behaviour 
has been recently demonstrated [16]. Formally, the figure 
of merit has no upper bound, so that arbitrarily large val- 
ues can be obtained within a model approach. Though 
for the realistic molecular junctions we are going to con- 
sider farther below the effective electrode-molecule cou- 
pling is on average strong (~100 meV), the combined use 
of semiconducting electrodes and chemical tuning still al- 
lows for thermopower optimization [ ]. 
Thermopower of single-molecule junctions— In this sec- 
tion we will discuss to which extent the analysis pre- 
sented in the previous paragraphs can be realized in re- 
alistic molecular junctions. Methodologically, we will 
use a density-functional parametrized tight-binding ap- 
proach (DFTB) [30, 31] combined with Green function 
techniques to calculate the transmission function of the 
junctions and from it the thermoelectric transport prop- 
erties at zero applied bias. We have addressed two main 
issues: (i) the role of the linker groups in determining 
the relative position of the HOMO with respect to the 
Fermi energy, and (ii) modifications of the thermopower 
through selective chemical doping of the molecules. 
Figs. 3 and 4 show the transmission spectra of several 
junctions together with the corresponding thermopower 
S and figures of merit ZT. From Fig. 3, where the influ- 
ence of different linker groups on the transmission func- 
tion around the Fermi level [ ] is shown, we see that thiol 
linkers turn out to be the most effective linkers. They in- 
duce a high-transmission HOMO resonance around 0.5 
eV below the Fermi level. Notice also the sharp sup- 
pression of the transmission tails above the conduction 
band edge due to the absence of spectral weight in the 
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FIG. 2: Top panel: Energy profile of the model 
semiconductor-molecule-semiconductor junction. Egap is the 
gap between the valence (VB) and conduction (GB) bands, 
A = ^HOMO — denotes the position of the HOMO with 
respect to ^f- The molecular gap is chosen so as to have the 
LUMO level placed inside the semiconducting gap. Lower 
panel: Seebeck coefficient S (b) and thermoelectric figure of 
merit (ZT) (c) as a function of A and the effective coupling 
strength to the electrodes F at T = 300 K. Notice that the 
largest thermopower could be achieved by having a weak cou- 
pling to the electrodes (large tunnel barriers) and a molecular 
resonance close to the Fermi level (large derivative). Since 
there is no upper bound for ZT, very large values can be 
formally attained within these model calculations. 



semiconducting electrodes: this effect can additionally 
lead to a large derivative at the Fermi level and to a 
dramatic increase in S when comparing with metallic 
electrodes. The other linkers studied here induce larger 
energy shifts of the transmission resonances away from 
the Fermi level and thus seem to be less appropriate for 
the realization of an efficient thermoelectric molecular 
junction. Taking this into account we have then investi- 
gated small polyacene molecules (benzene, napthtalene, 
anthracene) with NH2 functional groups added at the 
ortho-positions of the benzene rings, see Fig. 4. As re- 
cently shown [20] , the controlled torsional motion of side 
groups can help to tune the thermopower of a molec- 
ular junction; this would require however very specific 
mechanical manipulations at the single molecule level. 
In the present study we demonstrate a different alter- 
native route to potentially achieve high thermoelectric 
efficiency via pure chemical tuning of transmission reso- 
nances in combination with non-metallic electrodes. As 
it comes out from the analysis of our results, the side 
groups act as a chemical gate which shifts the transmis- 
sion spectrum towards the Fermi level. As a result, for 
the three polyacene species studied the HOMO levels are 



3 




150 250 350 

Temperature (K) 



50 150 250 3! 

Temperature (K) 



FIG. 3: Top panel: Zero-temperature transmission spectra 
of a benzene molecule sandwiched between Si electrodes via 
different linker groups. The position of the HOMO level is 
modified due to charge transfer effects. Bottom panel: Tem- 
perature dependence of the thermopower S (left) and corre- 
sponding figure of merit ZT (right) for the different linker 
groups. 
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FIG. 4: Top panel: Zero-temperature transmission spectra 
of three polyacene molecules (benzene, naphtalene, and an- 
thracene) with NH2 side groups and attached to Si electrodes 
via thiol linkers. For all three molecules the HOMO level 
lies almost at the corresponding Fermi energy. Notice also 
that for the benzene molecule the HOMO state moves into 
the semiconductor band gap. Bottom panel: Temperature 
dependence of the thermopower S (left) and figure of merit 
ZT (right). 



essentially lying at or slightly below Ep. Moreover, due 
to a low weight onto the linker region, the coupling of 
these levels to the electrodes is rather weak, leading to 
a strong narrowing. Notice that in contrast to naphtal- 
ene and anthracene, the HOMO resonance of the benzene 
dithiol molecule has a very low transmission; this is re- 
lated to the fact that the energy shift induced by the side 
groups moves this level into the semiconducting band 
gap, where no spectral support exist in the electrodes. 
Taking into account the narrowness of the HOMO levels 
together with their energetic position, we may expect a 
dramatic increase of the thermopower when comparing 
with the junctions studied e.g.^ in Fig. 3; this is the case 
as shown in the lower panel of Fig. 4, where maximum S- 
values ranging from 400—1000 fiV/K could be achieved. 
Also very large figures of merit of up to 60 at about 270 
K might be ideally reached. The fact that the maximum 
values of S and ZT do not necessarily lie at the same 
temperature demonstrates that the optimization of the 
thermopower alone is not enough to achieve optimal fig- 
ures of merit. We stress here that our main focus are 
electronic contributions while the vibrational component 
of the thermal conductance has been neglected at this 
stage. Experiments [13, 14] as well as theoretical stud- 



ies [ ] suggest that vibrational contributions in organic 
molecular systems are of the order of pW/K above room 
temperature. We expect that the high mismatch between 
the silicon surface modes and the local vibrational modes 
of the molecules will lead to similar orders of magnitude 
in our case. On the other hand, our calculated electronic 
thermal conductance for e.g., the naphtalene molecule 
is about 1 pW/K at 1000 K. Notice that we can write 
ZT = ZTei /^ei/(^ei + ^ph)- Hcucc, the inclusion of vi- 
brational contributions would lead to a reduction of the 
second factor by at most one order of magnitude. Still, 
the resulting ZT values (~ 4-6) would be large when 
compared to nanowires [ , ]. 

In summary, we have shown in this study that ther- 
mopower engineering in molecular junctions could be at- 
tained by an appropriate combination of chemical gating 
and the use of semiconducting electrodes. The latter pro- 
vide a way to eliminate the LUMO channel from the ther- 
mopower, thus lifting the partial cancelation of hole- and 
electron-like contributions. Though the obtained very 
large figures of merit should be considered with caution 
due to the neglect of vibrational contributions, large ZT 
values might still be achieved by using molecular junc- 
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tions. We also expect that our results can be scaled up 
to estimate the thermoelectric efficiency of self- assembled 
monolayers if the molecule-molecule interactions are not 
very strong. 
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